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Abstract 

We present a dynamical study of e + e~ production in C + C and Ca + Ca 
collisions at BEVAL AC /SIS energies on the basis of the covariant transport 
approach HSD employing momentum-dependent p-meson spectral functions 
that include the pion modifications in the nuclear medium as well as the 
polarization of the p- meson due to resonant p—N scattering. We find that 
the experimental data from the DLS collaboration cannot be described within 
the />meson spectral function approach. A dropping 77-mass scenario leads 
to a good reproduction of the DLS dilepton data, however, violates the mr- 
scaling of 7r° and r\ spectra as observed by the TAPS collaboration as well as 
r\ photoproduction on nuclei. 
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I. INTRODUCTION 



The mechanism of chiral symmetry restoration at high baryon density is of fundamental 
importance for the understanding of QCD |jJ0], but a clear experimental evidence has not 
been achieved, so far. In this respect, dileptons are particularly well suited for an inves- 
tigation of the early phases of high-energy heavy-ion collisions because they can leave the 
reaction volume essentially undistorted by final-state interactions. Indeed, dileptons from 
heavy-ion collisions have been measured by the DLS collaboration at the BEVALAC 
and by the CERES 0§, HELIOS 0JT§, NA38 ]ITJ and NA50 collaborations at SPS 
energies. 

Various theoretical calculations HT3HI7H fail in describing the spectral shape of the dilep- 
ton yield for nucleus-nucleus collisions at SPS energies when employing a 'free' p-meson 
form factor but are successful when using a dropping p-mass in the medium. However, a 
more conventional approach including the change of the p-meson spectral function in the 



medium due to the coupling of p, tt, A and nucleon dynamics along the lines of Refs. [|18| -|2T 



was found to be (roughly) compatible with the CERES data [0,^T|] as well, whereas the 
dimuon data of the HELIOS-3 collaboration || could only be described satisfactorily when 

Meanwhile, 



including 'dropping' meson masses p2| , p3| according to Hatsuda and Lee 
our knowledge on the p spectral function has improved since - as first pointed out by Friman 



and Pirner |25[ - resonant p-N interactions in p-wave scattering significantly enhance the 
strength in the vector-isovector channel at low invariant mass. In fact, the CERES data 
for S + Au at 200 A GeV and Pb + Au at 160 A GeV, using an expanding fireball model, 
were found to be compatible within such a hadronic scenario p6[ . This general shape of the 



p spectral function has also been confirmed recently in Ref. |27j where the authors find a 
significant smearing of the p strength at low p momenta due to a selfconsistent evaluation 
of the resonance widths. Furthermore, as shown in Refs. [p8,29[], both scenarios (dropping 
mass and p-meson broadening) are compatible with QCD sum rules. 

The thermodynamical result of Ref. p6[ was supported by HSD transport calculations 
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30| where an improved p spectral function {i.e. by constraints from 7p and 7 A data) was 
involved. The important conclusion in Ref. was that both the 'dropping mass' sce- 
n],[E|^,[n|] as well as the hadronic spectral function approach |26| lead to dilepton 



nario 



spectra which are in good agreement with the experimental data for all systems at SPS 
energies. 

In this respect, it is a natural step to investigate within the same framework dilepton 
production at BEVALAC energies where quite a different temperature and density regime 
is probed. From a recent publication of the DLS collaboration ||, it turns out that the new 
measurements of dilepton yields in A + A collisions have increased by about a factor of 6-7 
in comparison to the previous data set |3|-|5| due to an improvement of the DLS detector 
and the data analysis correcting for dead time losses. However, the first generation of DLS 



data were reasonably described by several theoretical models [|32| -|35|j without incorporating 
any medium effects. These various models are consistent in the general tendency that the 
dilepton pairs of invariant mass below 0.4 GeV are primarily produced from hadronic sources, 
such as the Dalitz decays of 7r°, A and rj and pn bremsstrahlung, whereas for M > 0.4 GeV 
the dominant contributions stem from the ir + n~ annihilation channel and direct decays of 
vector mesons (p, u). Since this has lead to a fair agreement with the old DLS data the 
aforementioned theoretical models are not able to explain the new DLS data set || without 
incorporating new ingredients such as vector meson selfenergies. 

The DLS collaboration pointed out that the inconsistency between the data and theoret- 
ical results - incorporating no medium effects - might be due to a strong underestimation of 
the r) contribution in the model calculations. As found in Ref. H the shape of the dilepton 



yield for M < 0.4 GeV is quite similar to that of the r\ component in the BUU model |35| 



The main knowledge about rj production in heavy- ion collisions at BEVALAC/SIS energies 
is provided by the TAPS collaboration p6|-|38[|. However, TAPS measures 7r° and rj mesons 
in a small rapidity interval near y cm whereas the DLS spectrometer covers the forward ra- 
pidity region. Within a thermal model calculation - using the TAPS mid-rapidity data and 
assuming an isotropic rj angular distribution - the DLS collaboration found that the rj yield 



in their rapidity window then is in good agreement with the BUU model predictions |35| , |39 . 



The dilepton enhancement observed by the DLS Collaboration thus might be due to a strong 
anisotropy of the rj angular distribution in the center-of-mass system. 

In this article we address the DLS 'puzzle' by employing in our transport calculation 
novel theoretical developments, i.e. in-medium p spectral functions as well as hypothetical 
in-medium modifications of the rj mass. To arrive at more stringent conclusions concerning 
the rj contribution we simultaneously analyze the DLS dilepton data || together with the new 
TAPS data with respect to ttit spectra of 7r°'s and 77's as well as r\ photoproduction 

on nuclei. 

Our paper is organized as follows: In Section 2 we briefly describe the transport approach 
employed as well as the elementary production channels. Section 3 contains a comparison 
of our calculations for 7r° and rj m^-spectra with the TAPS data. In Section 4 we discuss 
the p spectral functions which will be applied for calculating the dilepton spectra. Section 5 
contains a detailed comparison of the calculated dilepton spectra at BEVALAC energies 
with the DLS data ||. In Section 6 we investigate the consequences of a dropping 77-mass 
scenario for m^-spectra and dilepton production at BEVALAC/SIS energies as well as 77 
photoproduction on nuclei. The 77 angular anisotropy is investigated in Section 7 while we 
close with a summary and discussion of open problems in Section 8. 

II. TRANSPORT MODEL INGREDIENTS 

In this work we perform the theoretical analysis along the line of the HSD approach 
f4Cfl which is based on a coupled set of covariant transport equations for the phase-space 
distributions fh(x,p) of hadron h p0| , |41| 



i.e. 



= Yl /rf2d3d4...[G t G']i 2 ^34..5r(n + n 2 -n3-n4...) 

h2h 3 h4... 

x { A 3 (z> Ps)fh 4 (x,p 4 ) fh (x, p) f h2 (x, p 2 ) 
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- fh(x,p)fh 2 (x,P2)fh 3 (x,P3)fh i (x,p 4 )}... ■ (1) 

In Eq. ([1]) U^(x,p) and Ufi(x,p) denote the real part of the scalar and vector hadron selfen- 
ergies, respectively, while [G + G]i2^34...£r(n + n 2 — n 3 — n 4 . . .) is the 'transition rate' for the 
process 1 + 2^3 + 4 + .... In quantum many-body systems this transition rate is partly off- 
shell - as indicated by the index T of the 5-function - with width T « r~^ n where r co u is the 
average collision rate for the colliding hadrons. In our semi-classical transport approach the 
on-shell limit r — > is adopted which is found to describe reasonably well hadronic spectra 
in a wide dynamical regime. Note that the on-shell limit has been proven explicitly to be 
a valid approximation for intermediate energy heavy-ion collisions in Ref. ||42|| . The hadron 
quasi-particle properties in (^) are defined via the mass-shell constraint f4~If , 5(ILjIF — M^ 2 ), 
with effective masses and momenta given by 

M* h (x,p) = M h + U*(x,p) 

W{x,p)=pl i -U£(x,p) , (2) 

while the phase-space factors 

f h (x,p) = 1 ±f h (x,p) (3) 

account for fermion Pauli-blocking or Bose enhancement, depending on the type of hadron 
in the final and initial channel. The dots in Eq. (TJ) indicate further contributions to the 
collision term with more than two hadrons in the final/initial channels. The transport 
approach ([J) is fully specified by U^(x,p) and U£(x,p) (p = 0,1,2,3), which determine 
the mean-field propagation of the hadrons, and by the transition rates G^G5^(. . .) in the 
collision term that describe the scattering and hadron production/absorption rates. For 
more details we refer the reader to Ref. [fHJ . 



In addition to the direct production of all mesons in Ref. |4(J we treat the production 
of vector mesons (p, u, (ft) and 77's at low energies also perturbatively. The vector mesons as 
well as 77 mesons are produced in baryon-baryon and meson-baryon collisions according to 
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the elementary production cross sections evaluated in Refs. [0,0 on the basis of boson 



exchange diagrams. For an explicit description of the parametrizations employed we refer 
the reader to Eqs. (14) and (18) and Tables 1-3 in Ref. [(43| for p, u and (j> production in 
pion-baryon and baryon-baryon channels. The rj cross section from pp collisions is taken in 
line with the boson exchange calculation in Ref. [11], which can be fitted as 



o-w-WVS) - a b + y Q- _ V ^L )2 (4) 
with a = 0.4 mb and b= 0.552 GeV 2 while y/so is the rj production threshold. This 



parametrization also reproduces the rj production data from the PINOT Collaboration f45 
rather well. The pn — > pnrj cross section is taken to be about 6 times the pp — > pprj cross 
section close to threshold in line with the new data from the WASA collaboration j|6) . The 



production of rj mesons in pion-nucleon collisions as well as rj scattering and reabsorption 
on nucleons is taken from Ref. [|7j employing the experimental data from [0] and using de- 



tailed balance as well as the known branching ratios of the iV(1535) to pions and 77-mesons. 
The respective rj absorption cross section was found in Ref. |1TJ (Fig. 2) to be well in line 
with the experimental rj photoproduction data on nuclei (see below). 

Apart from the production cross sections of the vector mesons p and u their final state 
interactions with nucleons have to be considered. For the p we adopt the resonance coupling 
model specified in Section 2 of Ref. [Q whereas for the ui we use the total and inelastic cross 
sections (Eqs. (34), (35)) of Ref. |51] which are less constraint by experimental data. The 



latter uncertainty, however, plays no essential role in the present analysis because the u- 
mesons dominantly decay in the vacuum and thus do not contribute to the dilepton invariant 
mass spectrum below 0.6 GeV. 

As noted above, the transport approach (HSD) was found to describe reasonably 
well hadronic data at SIS energies as well as hadronic and dilepton data at SPS ener- 



gies imramisD 
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III. M T SCALING 



In relativistic heavy- ion collisions at BEVALAC/SIS energies the nuclei can be com- 
pressed up to about 3 times normal nuclear matter density po- I n the hot compression zone 
the nucleons are exited to baryonic resonances which decay by emitting mesons and the 
produced mesons then can be absorbed, re-emitted and re-scattered. 

Information about the degree of equilibration can be obtained from itlt spectra, where 
rriT = (Pt + m 2 ) 1 / 2 and px is the transverse momentum of a particle with mass m. In this 
case the spectrum 1 j 'm\da / 'dm? is of Boltzmann type ~ exp(— /3m^) with a slope parameter 
(5 being related to the global (inverse) temperature at freeze-out in the absence of flow. 

The transverse-mass spectra of ir° and i] mesons in heavy-ion collisions at SIS energies 
were recently measured by the TAPS collaboration p6|-|38|| where a mx scaling has been 
observed for both mesons. Such a universal property of the meson spectra at SIS energies 



has been predicted by the Quark Gluon String Model calculations in Ref. |52| for the Ar + Ca 
system. In the remainder of this section we will thus systematically investigate itit spectra 
for nucleus-nucleus collisions at SIS energies within the transport approach specified above. 

In Fig. [I] we compare the results of our calculation for the transverse-mass spectra of 
7T° and rj mesons with the TAPS data. The upper part shows the my spectra for 7r°'s 
(dashed histogram) and 7/'s (solid histogram) for C + C at 1.0 A-GeV in the rapidity interval 
0.42 < y < 0.74 and at 2.0 A-GeV for 0.8 < y < 1.08. The experimental data - the open 
circles and solid squares correspond to tt° and 77 mesons, respectively - are taken from 
Ref. [p7H . The theoretical results as well as the experimental data at 2.0 A-GeV here are 
multiplied by a factor of 10 2 . The middle part corresponds to Ca + Ca at 1.0 A-GeV for 
0.48 < y < 0.88 (multiplied by 10 _1 ) and at 2.0 A-GeV for 0.8 < y < 1.1 in comparison 
with the data from Ref. ||38|| . The lower part shows the calculated spectra for Ni + Ni 
at 1.93 A-GeV for 0.8 < y < 1.1 in comparison with the data from Ref. [38[ . 



As seen from Fig. [I], the HSD transport model gives a reasonable description of the itit 
spectra of pions and 77's as measured by the TAPS collaboration without incorporating any 
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medium modifications for both mesons. It is important to point out that these calculations 
are parameter-free in the sense that all production cross sections for r\ mesons are extracted 
from experimental data in the vacuum and the 7/-nucleon elastic and inelastic cross sections 
are obtained by using detailed balance on the basis of an intermediate N(1535) resonance. 



IV. IN-MEDIUM p PROPAGATOR 

The dilepton production is calculated perturbatively by including the contributions from 
the Dalitz-decays A -> Nl + l~, N* -> Nl+l~, vr° -> 7/+/-, r, -> , u -> n°l+r and 

the direct dilepton decays of the vector mesons p, uo and 0, where the p and meson may 
be produced in 7r + 7r~ and KK or up collisions, respectively. For a detailed description of 
A, N* Dalitz decays we refer the reader to Ref. ||33|| ; the dilepton decays of the rj and vector 



mesons (u,<p) are described in Ref. 31 



The dilepton radiation from p mesons is calculated as 

dN l+l - 2M 

——— = -Br lmD p (q , q; p B , T), (5) 

CLM 7T 

where D p is the p-meson propagator in the hadronic medium depending on the baryon 
density ps and temperature T as well as on energy go and 3-momentum q = |gj in the 
local rest frame of the baryon current ('comoving' frame) while Br ~ 4.5 x 10~ 5 is the 
experimental branching ratio to dileptons. The invariant mass M is related to the p-meson 
4-momentum in the nuclear medium as 

M 2 = ql - g 2 , (6) 

while 

\mD p = l -{\mD L p +2\mD T p ) (7) 

is spin averaged over the longitudinal and transverse part of the p-propagator, 
t nL,T, rp\ ImS^ T (g 0) q- p B , T) 

ImD (g , g; Pb, T) = p 2 — . (8) 

| M 2_( m (0)^ _ E ^( g0)?;pB)T )|2 
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Following Refs. p6[ , |5"3f the p-meson selfenergy Y p (qo,q, Pb,T) is obtained by combining 



the effects of the different hadronic interactions: 



Y p ' T — E^ + ^ pBB -i + + ^pKK-L- (9) 



The explicit evaluation of the various selfenergy components is discussed in detail in Ref. p6 
We here employ an improved version in the following respects: 

(a) in addition to p-wave pN — > B interactions (B=N, A, iV(1720) and A(1905)) we also 
account for s-wave excitations into iV(1520), A(1620) and A(1700) resonances (the 
corresponding coupling constants usual, estimated from the pN partial decay 
width) ; 

(b) the coupling of the resonances to (virtual) photons is calculated within the improved 



vector dominance model of Kroll, Lee and Zumino 541 to avoid overestimates of 



the B —>■ Nj branching ratios. This leads to a modification in the coupling of the 



transverse part of the p propagator entering Eqs. (|5|), (|TT|) such that the combination 



(mjM) 4 i m £)J(gr 0) q- p B) T) is replaced by the 'transition form factor' 

F T (q , q; p B ) = -Im[E^ + E^ ai + Y? pKKl ] \d p - 1| 2 - ImSj BB _i \d p - r B \ 2 

M 2 - Y T -Y T -Y T -r n Y T 
1VI pirn ^p-Ka x l -'pKK 1 'B^pBB- 1 

M 2 - (m P 0) ) 2 - Ej 



d M^ q; pb) = pn \^ pna ; {oh r K L " > (io) 



where 

Pb 



fpBN \rn> p 'Y 



denotes the ratio of the photon coupling to its value in the naive VDM. 
(c) the medium modifications of the two-pion selfenergy E p7r7r are extended to arbitrary 



3- momentum within the recently developed model of Urban et al. p5] 



The combined p-meson selfenergy is then further constrained by experimental information 
on 7p and 'jA absorption cross sections as described in Ref. ||53|| . With the aforementioned 
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improvements (a) to (c) a satisfactory description of the photoabsorption data, which rep- 
resent the M 2 — ► limit of the dilepton regime, can be achieved on protons as well as on 
nuclei. In fact, the melting of the resonance structures (above the A mass), as seen in the 



photoabsorption data on nuclei ||56|| , can be explained by the broadening of the higher reso- 
nances due to a strong coupling to the short-lived p-meson in the medium ]53 . This provides 
a proper basis to assess the dilepton radiation originating from in-medium nn annihilation 
and p decays in A + A collisions at BEVALAC/SIS energies. 

As an example - relevant for the heavy-ion reactions at BEVALAC/SIS energies - we 
show in Fig. || the spin averaged — ImD p (g , q\ Pb, T) in the approach of Ref. |26j (including 
the modifications (a) to (c)) as a function of the invariant mass M and the 3-momentum 
q for a temperature of 70 MeV at p# = 0,p ,2p an d 3p , respectively (note that at this 
temperature about 10% of the nucleons are excited into A's such that p^ — 0.9p#). With 
increasing baryon density we find the p spectral function to increase substantially in width 
showing only minor structures at high density. Thus the lifetime of the p-meson in the 
nuclear medium becomes very short due to well established hadronic interactions. 

Furthermore, we will also employ the p spectral function from Ref. |27| which is calculated 
by taking only the resonance-holes states (without the correction Eq. (|TUD) into account 
including, however, all nucleon resonances up to a mass of ~ 1.9 GeV. In addition, in the 
latter calculation the feedback of the in-medium p spectral function on the resonance widths 
is taken into account selfconsistently, i.e. iterated out at given nuclear matter density and 
3-momentum. In Fig. [3] we show for comparison the spectral function from Ref. |27j for 
p = 0, po)2po an d 3po at temperature T = 0. While the overall structure is similar to that 
in Fig. [| the spectral functions in the approach of Ref. [26] are smoother at low invariant 
mass and the structures at small momenta q, which result from a selfconsistent treatment 



of the resonance widths in Ref. ||27|| , are washed out. 

The cross section for the pion annihilation channel tt + tt~ — > p° — > e + e~ is taken in line 
with Ref. i57H 



as 
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<W_ef e- (M) = ^ K 0) ) 4 ^D p (q , g; p B) T) (11) 

where fc = (M 2 — 4m 2 r ) 1/ ' 2 /2 is the pion 3-momentum in the center-of-mass frame and a is 
the fine structure constant. The pmr coupling constant and bare p mass mffl are fixed to 
reproduce p-wave nir scattering and the pion electromagnetic form factor in free space |]26| , |27 . 



Our dynamical calculations are carried out as follows: the time-evolution of the proton- 
nucleus or nucleus-nucleus collision is described within the covariant transport approach 



HSD |4yj31|] without any dropping vector meson masses. Whenever a p-meson is produced in 
the course of the hadronic cascade (by baryon-baryon, meson-baryon or pion-pion collisions), 
its 4-momentum in the local rest frame of the baryon current is recorded together with the 
information on the local baryon density, the local 'transverse' temperature and its production 
source. We note that the definition of a local temperature is model dependent; here we have 
used a logarithmic fit to the transverse pt spectra of mesons at midrapidity. Without going 
into a detailed discussion of this issue we note that our results for dilepton spectra do not 
change within the numerical accuracy when using a constant temperature T = 70 MeV at 
BEVALAC energies since lmD p depends rather weakly on T (at fixed nucleon density). 



V. DILEPTON SPECTRA IN COMPARISON TO BEVALAC DATA 

In Fig. |]we present the calculated dilepton invariant mass spectra da/dM for Ca + Ca 
(upper part) and C + C (lower part) at the bombarding energy of 1.0 A-GeV integrated 
over impact parameter and compare them with the new experimental data of the DLS 
collaboration || including the DLS acceptance filter (version 4.1) as well as a mass resolution 
AM/M = 10%. The thin lines indicate the individual contributions from the different 
production channels; i.e. starting from low M: Dalitz decay tt° — ► 7e + e~ (dashed line), 
r\ — ► 7e + e _ (dotted line), A — * Ne + e~ (dashed line), uj — > 7r°e + e~ (dot-dashed line), iV* — > 
Ne + e~ (dotted line), proton- neutron bremsstrahlung (dot-dashed line), ttN bremsstrahlung 
(dot-dot-dashed line); for M ~ 0.8 GeV: uj e + e~ (dot-dashed line), p° — > e + e~ (dashed 
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line), 7r + 7r~ — > p — > e + e~ (dot-dashed line). The pion annihilation channel as well as 
the direct decay of the vector mesons here have been calculated with the 'free' p spectral 
function. The full solid line represents the sum of all sources. 

Our present result is consistent with that from Ref. |35[ at low invariant mass M < 
0.4 GeV. However, the contributions from direct decays of p and u mesons are somewhat 
larger because novel elementary production cross sections for pion-baryon collisions from 
Ref. [^H were taken into account as well as p and u elastic and inelastic scattering processes 



with nucleons using the cross sections from Ref. |p0| , pl| . As seen from Fig. (| the new 
BEVALAC data cannot be properly described in terms of a 'free' p spectral function. The 
discrepancy between the data and the calculations for Ca+Ca as well as for C+C at 0.15 < 
M < 0.4 GeV is about a factor of 3-5. At M ~ m p the calculation is within the error bars 
except for the last experimental point. 

The description of the DLS data is slightly improved when including the p spectral 
functions from Ref. [^6j (cf. Fig. |^) or Ref. [[Z7J (cf. Fig. |3|). In Fig. ^ we show the dilepton 
spectra (full solid line) for Ca + Ca (upper part) and C + C (lower part) at 1.0 A-GeV 
calculated with the p spectral function from Ref. (Fig. 0) in comparison with the data || . 
The same comparison is carried out in Fig. |^ for the spectral function (at T=0) from Ref. [27 



(Fig. |H). The dashed lines correspond to the channel p° — > e + e~ while the dot-dashed lines 
indicate the contribution from 7r + 7r~ — > p — ► e + e~ as described in Section 4. The full p 
spectral functions lead to a shift of the pion annihilation contribution to the small invariant 
mass region as well as to a broadening of the p-meson contribution in both approaches. 
We note that in our transport approach the pions are on-shell, thus we can not probe the 
p-spectral function in the reactions 7r + 7r~ — > p° — > e + e~ below the two pion threshold at 
M < 2m 7r . The tails from the -k + tt~ annihilation and p meson at M < 2m n , which are seen 
in Figs. [| and |6|, exist only due to the finite mass resolution of AM/M ~ 10%. We note 
that we discard meson-baryon (mB) and meson-meson (mm) bremsstrahlung channels as 
well as the Dalitz decays of the baryon resonances (stemming from secondary pion induced 
reactions) in order to avoid double counting when employing the full p spectral functions. 
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Due to the on-shell treatment of pions in the transport approach the dilepton spectra for 
M ~ 2m v might not be properly described. In order to investigate the effects from an off- 
shell propagation of pions in the channel tt + tt~ — ► p° — > / + /~, we have performed additional 
calculations within the thermodynamical approach used in Ref. [^6] at SPS energies. For 
BEVALAC energies we use exactly the same procedure, however, with temperature and 
density profiles obtained from the HSD transport model (|0| for Ca + Ca at 1.0 A-GeV. 
The solid curves in Fig. [7] correspond to the sum of all channels and the p decay (from 
7r + 7r~ annihilation and direct production), respectively, while the dashed lines are from the 
transport model (cf. Fig. [5]). As seen from Fig. [5] an off-shell pion propagation leads to 
a large enhancement in the p-decay channel especially below the 27r threshold. However, 
considering all channels simultaneously, the maximum increase - compared to the transport 
model - is only 30%. Thus, all independent calculations give practically the same results; 
the new DLS data are underestimated at least by a factor of 3 for 0.15 < M < 0.4 GeV even 
when taking into account the modification of the p meson properties in the nuclear medium. 



VI. DROPPING i] MASS? 

According to the general scaling idea of Brown and Rho 77 mesons might also change 
their properties in the medium. In this section we examine the possibility of a dropping 77 
mass at BEVALAC/SIS energies as well as for rj photoproduction on nuclei. For our analysis 
we use a linear extrapolation of the rj mass, 

m; = m° v (l-a^y (12) 

with a 'large' value for a m 0.18, which might be considered as an upper limit. According 
to Eq. fll2|) the 77 mass decreases at normal nuclear matter density by 18%. The maximum 
baryon density reached in central Ca + Ca collisions is about 2.5po, thus the dropping 77 
mass leads to an essential reduction of the 77 production threshold in BB and mB collisions 
and to a strong enhancement of the 77 population in heavy-ion collisions. Similar, but less 
pronounced effects are expected for the 77 photoproduction on nuclei for densities ps < Po- 
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A. Nucleus-nucleus collisions 

In Fig. |]we show the dilepton spectra (full solid line) for Ca + Ca (upper part) and C + C 
(lower part) at 1.0 A-GeV calculated with the full p spectral function |26| and a dropping 77 
mass in comparison with the data from Ref. II. As seen from Fig. [8], the dropping 77 mass 
scenario together with the full spectral function approach leads to a good reproduction of 
the BEVALAC data. 

However, the dropping rj mass scenario does not yield the observed ttit scaling. The 
transverse-mass spectra of 7r° and 77 mesons for C + C in 0.42 < y < 0.74 (lower part) and 
Ca + Ca for 0.48 < y < 0.88 (upper part) at 1.0 A-GeV are shown in Fig. |. The assignment 
of lines is the same as in Fig. [l]. The dot-dashed histograms correspond to the 77 m^-spectra 
obtained within the dropping r\ mass scenario. The rnx scaling in this case is violated by 
about a factor of 3 especially at low my ~ m v . 

Thus we find that the simple 77 dropping mass scheme is not consistent with the my 
scaling observed by the TAPS collaboration. A more sophisticated scenario for the 77 self- 
energy (like momentum dependent 77 potentials etc. |59[]) still has to be studied. 

B. i] photoproduction on nuclei 

Since the dynamics of nucleus-nucleus collisions at SIS energies is rather complex it is 
useful to have further constraints on the 77 properties in the medium as e.g. tested in 77 
photoproduction on nuclei where the kinematical conditions are much better under control. 
The latter problem has been addressed in detail in Refs. ||],fnj on the basis of transport 
approaches, where further details of the explicit calculations can be found. Since in the 
present study we want to explore the possible consequences of a dropping 77 mass at finite 
nuclear density we reparameterize the 77 photoproduction on a proton as 

v-yp^rip(Vs) = a(l - x) b exp - v^o) 2 / c ) ( 13 ) 
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with x = yfso/^,a = 161.5 fib,b = 0.571 and c = 8.9 10 3 GeV 2 , which provides a 



good description of the experimental data from Refs. p2lp3|. The rj photoproduction on a 



neutron is found experimentally to be 2/3 as that on a proton; we adopt the same ratio 



for our calculations. The parametrization (O) as a function of the invariant energy above 



threshold a/s — ^/sq now allows to simulate r] photoproduction with a dropping rj mass 
according to Eq. (|i~2"[) . We note, that the parametrization ([0|) assumes the photoproduction 
to be dominated by phase space or the invariant energy above threshold; this assumption 
still has to be controlled by microscopic in-medium calculations on rj photoproduction and 
here only serves as an estimate. 

The results of our transport calculations for rj production on 40 Ca for photon energies 
from 600 - 780 MeV are displayed in Fig. [10] in comparison to the data from Ref. [p9j] . The 



full line corresponds to the bare 7] mass whereas the dashed line is obtained for a dropping 
mass according to Eq. (|T2|) with a = 0.18. The excitation function clearly indicates that the 
r) properties do not change much in the medium except for a strong rj absorption which is 
about 65 % in case of 40 Ca (cf. Ref. |61|); the dropping //-mass scenario leads to enhanced r\ 
cross sections especially at low photon energy which are not observed experimentally. Thus 
also rj photoproduction data exclude a dramatic change of the r/ properties in the medium 
as possibly 'requested' by the DLS data (cf. Fig. |]). 

VII. ANISOTROPY OF THE rj ANGULAR DISTRIBUTION 

As proposed by the DLS collaboration || the inconsistency between the DLS and TAPS 
data might be due to differences in the experimental acceptance - mid-rapidity for TAPS 
versus forward rapidity for DLS - and a strong anisotropy of the r\ angular distribution. The 
TAPS group calculated the total 7] cross section from a model where mesons are emitted 
isotropically from a thermal source at y cm . Using the mid-rapidity TAPS data and a simple 
thermodynamical model the DLS collaboration reproduces the result from our previous 
calculation . 
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In Fig. |TT| we, therefore, present the result of our transport calculation for the angular 
distribution of rj mesons for Ca + Ca at 1.0 A-GeV. Mainly due to rescattering in matter 
the rj angular distribution is anisotropic; however, this anisotropy is not very pronounced. 
As was estimated by the TAPS collaboration, the implementation of this anisotropy gives 
at maximum a 20% correction for the total rj cross section [ |B3{] . Thus, the anisotropy effect 
according to our transport calculations cannot explain an enhancement by a factor of 3-7 in 
the new DLS dilepton data. 

VIII. SUMMARY 

On the basis of the covariant transport approach HSD ^(| we have studied dilepton 
production and rnj- scaling in C + C and Ca + Ca collisions at BEVALAC energies. We 
have found that the my scaling of 7r° and rj observed by the TAPS collaboration is reproduced 
by the transport calculation even for the small system C + C when incorporating no medium 
effects for the pions and etas. 

Various contributions are taken into account for dilepton production: the Dalitz-decays of 
A, A* resonances and ir°,ri,uj mesons, proton-neutron bremsstrahlung, nN bremsstrahlung 
as well as the direct dilepton decays of the vector mesons p and uo. It was shown that the 
new DLS data are underestimated by a factor of 6-7 in the transport calculations for 
invariant mass 0.15 < M < 0.4 GeV when no medium effects are involved as in our previous 



calculations 35 



Including the p meson modifications in the medium (according to the spectral functions 



from Refs. p6| , p7[ ) slightly improves the agreement with the data, but the discrepancy is 
still about a factor of 3. This also holds for the thermodynamical approach |26| where 
an off-shell propagation of pions is taken into account. A simple dropping rj mass scheme 
together with the full p spectral function can lead to dilepton spectra that are in a good 
agreement with the DLS data. However, the my-scaling of 7r°'s and 77's observed by the 
TAPS collaboration cannot be reproduced within this dropping rj mass scenario anymore. 
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The dropping mass scenario also is clearly incompatible with the r/ photoproduction data on 
nuclei. Also the anisotropy of the rj angular distribution due to rescattering effects in matter 
- as obtained from our transport calculation - is much too small to explain the enhancement 
of the dilepton yield at 0.15 < M < 0.4 GeV as observed by the DLS collaboration. 
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FIGURES 

FIG. 1. The calculated transverse-mass spectra of 7r° and r\ mesons in comparison with the 
TAPS data. The upper part shows the mx spectra for 7r°'s (dashed histogram) and r/'s (solid 
histogram) for C + C at 1.0 A-GeV in the rapidity interval 0.42 < y < 0.74 and at 2.0 A-GeV for 
0.8 < y < 1.08. The experimental data - open circles and solid squares correspond to 7r° and r\ 
mesons, respectively - are taken from Ref. The theoretical results as well as the experimental 
data at 2.0 A-GeV are multiplied by a factor of 10 2 . The middle part corresponds to Ca + Ca 
at 1.0 A-GeV for 0.48 < y < 0.88 (multiplied by 10" 1 ) and at 2.0 A-GeV for 0.8 < y < 1.1 in 
comparison with the data from Ref. [p£}|. The lower part shows the calculated tut spectra for 



Ni + Ni at 1.93 A-GeV at rapidities 0.8 < y < 1.1 in comparison with the data from Ref. [ 38 1 



FIG. 2. The (negative) imaginary part of the p propagator averaged over the longitudinal 
and transverse components (where the latter contain the correction Eq. ([H])) as a function of 
the invariant mass M and the momentum q for baryon densities of 0, 1 po, 2 po, and 3 po and 
temperature T = 70 MeV in the extended approach of Ref. [26]. Note the different absolute scales 
in the individual figures. 

FIG. 3. The (negative) imaginary part of the p propagator averaged over the longitudinal and 
transverse components as a function of the invariant mass M and the momentum q for baryon 
densities of 0, 1 po, 2 po, and 3 po at temperature T = MeV in the approach of Ref. j27j. Note 
the different absolute scales in the individual figures. 
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FIG. 4. The dilepton spectra (full solid line) for Ca + Ca (upper part) and C + C (lower part) at 
1.0 A-GeV calculated with the 'free' p spectral function including the DLS acceptance filter (version 
4.1) as well a mass resolution AM/M = 10% in comparison with the data from Ref. f§. The thin 
lines indicate the individual contributions from the different production channels including the 
DLS acceptance and mass resolution; i.e. starting from low M: Dalitz decay 7r° — > ^e + e~ (dashed 
line), rj — ► je + e~ (dotted line), A — > Ne + e~ (dashed line), ui — ► 7r°e + e~ (dot-dashed line), 
N* — * Ne + e~ (dotted line), proton-neutron bremsstrahlung (dot-dashed line), ttN bremsstrahlung 
(dot-dot-dashed line); for M « 0.8 GeV: u> — > e + e~ (dot-dashed line), p — ► e + e~ (dashed line), 
7r + 7r~ — > p — > e + e~ (dot-dashed line). 

FIG. 5. The dilepton spectra (full solid line) for Ca + Ca (upper part) and C + C (lower part) 
at 1.0 A-GeV calculated with the full p spectral function in the extended approach of Ref. |2(| 
in comparison with the data from Ref. || . The assignment of the individual contributions is the 
same as in Fig. |j. 

FIG. 6. The dilepton spectra (full solid line) for Ca + Ca (upper part) and C + C (lower part) 



at 1.0 A-GeV calculated with the full p spectral function in the extended approach of Ref. [27] 
in comparison with the data from Ref. || . The assignment of the individual contributions is the 
same as in Fig. |j. 

FIG. 7. The dilepton spectra for Ca + Ca at 1.0 A-GeV calculated within the thermodynamical 
approach |26| (solid lines) and within the transport model (dashed lines). 

FIG. 8. The dilepton spectra (full solid line) for Ca + Ca (upper part) and C + C (lower part) 
at 1.0 A-GeV calculated with the full p spectral function in the extended approach of Ref. |26| for 
the dropping r\ mass scenario in comparison with the data from Ref. Q. The assignment of the 
individual contributions is the same as in Fig. ||[ 
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FIG. 9. The transverse-mass spectra of 7r° and 77 mesons for Ca + Ca in the rapidity interval 
0.48 < y < 0.88 (upper part) and C + C for 0.42 < y < 0.74 (lower part) at 1.0 A-GeV. The 
assignment is the same as in Fig. [l]. The dot-dashed histograms correspond to the r\ m^-spectra 
obtained within the dropping rj mass scenario. 

FIG. 10. The calculated t] photoproduction cross section on 40 Ca as a function of the photon 



energy in comparison to the data from [49|. The solid histogram displays our results for a bare 



r\ mass whereas the dashed histogram corresponds to the dropping rj mass scenario fll2|). 

FIG. 11. The calculated angular distribution of rj mesons for Ca + Ca at 1.0 A-GeV in the 
center-of-mass system. 
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